Study of the electron-phonon interaction in metal diborides MeB2 (Me=Zr, Nb, Ta, 

Mg) by point-contact spectroscopy 
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We review investigations of the electron-phonon interaction (EPI) in metal diborides MeB2 
(Me=Zr, Nb, Ta, Mg) by point-contact (PC) spectroscopy. For transition metal compounds the 
PC EPI functions were recovered and EPI parameter A ^ 0.1 were estimated. The data are consis- 
tent with the measured surface phonon dispersion curves. The low A value questions some reports 
about superconductivity in these compounds. Contrary, EPI in superconducting MgB2 films man- 
ifests also in the PC spectra itself by virtue of an elastic EPI contribution to the excess current 
determined by the energy dependence of the superconducting order parameter. To analyse the 
phonon features in the PC spectra of MgB2 a two-band model is exploited and the proximity effect 
in the k-space is suggested. 

PACS numbers: 63.20.Kr, 74.80.Fp, 73.40.Jn 

Introduction. Recentlw discovered superconductivity in 
the sp compound MgB2tJ attracts much attention mainly 
due to relatively high Tc ~ 39 K, which is the highest 
as for two-component systems. Probably a specific fea- 
ture of this compound appears to be a rare example of 
two disconnected bands of th&.Eermi surface with quite 
different dynamical propertiead'Elu. One of those bands 
is two-dimensional (2D), with extremely strong electron- 
phonon iiji±|ejca|etion (EPI), while the other is 3D with a 
weak EPmeiy. 

The overwhelming majority of the communitV|| 
ipated an electron-phonon mechanism in MgB2Qi 
Thus, the preferable reason of the formation of Cooper 
pairs is EPI, and its investigation to understand the na- 
ture and peculiarities of superconductivity in MgB2 is 
vastly desirable. 

In the wake of the finding of superconductivity in 
MgB2 scientific activity was applied to its search in other 
diborides (all having a hexagonal layered crystal struc- 
ture of an AIB2 type) along with the further study of 
their properties. According to a recent reviewo no su- 
perconducting transition has been observed so far in the 
diborides of transition metals MeB2 (Me=Ti, Zr, Hf, V, 
Cr, Mo). As to the superconductivity in ]ii|tiB2, TaB2 
and ZrB2 controversial reports can be foundollj. 

In this article we would like to summarize our find- 
ings in the determination of the EPI function for the 
above mentioned diborides from the point-contact (PC) 
spectra. The attractive feature of the PC spectroscopy is 
that the measurement at low temperatures of the voltage 
derivative of the PC resistance allows a straightforward 
determinationB of the PC EPI function a^F{ijj). The 
knowledge of a^F{uj) for conducting systems provides a 
consistent check for the possibility of a phonon-mediated 
pairing mechanism, e.g., by estimation of the electron- 



phonon-coupling strength characterized by the EPI pa- 
rameter X = 2 J a'^F{uj)u!~^du). From the comparison 
of the experimentally determined a'^F{LLj) with theoreti- 
cal calculations different models and approaches can be 
discriminated. 

Experimental details. The samples are c-axis oriented 
filma-jof MgB2, whose characteristics are described in 
Ref.y and single crystals of MeB2 (Me=Zr, Nb, Ta). The 
residual resistivity po and the residual resistivity ratio 
(RRR) of all investigated diborides are shown in Table 1 . 

Different PC's were established in situ by touching of 
the MgB2 film or the cleaved surface of a transition metal 
single crystal by a sharpened edge of an Ag counterelec- 
trode. This geometry corresponds to a current fiowing 
preferably along the c-axis of the film. A number of con- 
tacts was measured by touching the film edge after break- 
ing AI2O3 substrate. By this means, the current flows 
preferably along the ab plane. As to PC's on another di- 
borides, the orientation of electrodes was not controlled 
with respect to the crystallographic axis. The zero- bias 
resistance Rq of the investigated contacts ranged from a 
few ohms up to several tens of ohms at 4.2 K. 

Both the differential resistance dV/dl and the second 
derivative oi I —V characteristic d'^V /dI'^{V) vs V were 
registered using a standard lock-in technique. 

Results and discussion. Accorjding to the Kulik, 
Omelyanchouk and Shekhter theorytj the second deriva- 
tive -d^I/dV^iV) of the / - curve of a balhstic 
PC at low temperatures due to EPI is proportional to 



^PC 



F{lj). In the free electron approximations 
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where R = dV/dl. The factor K = 1/2(1 - 61/ tan 6^), 
entering in apQ, takes into account the kinematic restric- 
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TABLE I: Parameters of investigated diborides. 



Compound 


Sample po 


■ fO^ flm 


R,R,R, 


n^-fO-^^, m-3 


MgBa 


film 


60 


2.3 


6.9 


ZrBa 


sgl. cry. 


3.3 


24 


f3 


NbBa 


sgl. cry. 






18.2 


TaBa 


sgl. cry. 


220 


f.2 


18.6 



" The density of carriers n was estimated by the number of 
valence electrons (2 for MgBa, 4 for ZrBa, 5 for NbB2 and 
TaBa) per volume of the corresponding unit cell. 



tion of electron scattering processes in PC (for transport 
and Eliashberg EPI functions the corresponding factors 
are: K = {1 — cos 6) and K—1, respectively, with 6 is 
the angle between initial and final momenta of scattered 
electrons). Respectively, the large angle 6 ^ n scatter- 
ing (back-scattering) processes of electrons dominate in 

According to Eq.(|lD the measured rms signal of the 
first Vi and second V2 harmonics of a small alternating 
voltage superimposed on the ramped dc voltage V defines 
the EPI function a|c(e) F{e): 



2 / 

(y-pc\ 



F(e) 



3V2 hvY V2 
~ edVf 



(2) 



To estimate the PC diameter d, which enters in Eqs.(|^) 
and (HLthe relation Rpc{T) ~ ^ + ^ derived by 
WexlertJ is commonly used, which consists of a sum of 
the ballistic, Sharvin, and the diffusive. Maxwell, terms. 
Here pi = pp/ne'^, where pp is the Fermi momentum, n 
is the density of charge carriers. 



A. Transition metal diborides 

Representative examples of the d^I/dV^{V) depen- 
dencies for each compound averaged over both voltage 
polarities are shown in Fig.|l| (left panel). Reproducible 
phonon maxima arc clearly resolved up to 100 mV (see 
ZrB2), while for TaB2 only low energy peak at 20 mV is 
seen. The spectra exhibit also a zero-bias anomaly, espe- 
cially remarkable in NbB2 and TaB2. A common feature 
for all spectra is the presence of the main low-energy max- 
imum placed at about 30, 28 and 20 mV for ZrB2, NbB2 
and TaB2, respectively. This is in line with the general 
consideration that at fixed spring constants the phonon 
frequency decreases with increasing of atomic mass (see 
Fig.|l], right panel). Such a behavior suggests that the 
first peak corresponds to the vibration of the transition 
metal. On the other hand the neutron data peak posL. 
tion for MgB2 on Fig.^ (right panel) at about 35meVE£l 
is far below the straight line connecting the MeB2 com- 
pounds. This might be considered as the consequence 
of a softening of the corresponding spring constants (i. 
e. metallic bonds in MgB2 instead of relatively strong 
MedB2p covalent bonds in the MeB2 series). 
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FIG. 1: Left panel: PC spectra -d^I(y)/dV^ oc 
d^V/dI^{dV/dI)~^ for investigated compounds at T = 4.2K. 
The dashed curves represent background behavior. The zero- 
bias resistance and modulation signal are for ZrBa: Ro = 
5.50, Vi(0)=0.8mV; for NbBa: Ro = 50 fi, "l/i(0)=2.8 mV; 
for TaBa: Ro = 25^1, Vi(0) = 1.3mV. Right panel: The po- 
sition of the first peak (squares) in the measured spectra for 
ZrB2, NbB2 and TaBa vs the inverse square root of atomic 
mass of the transition metals. For MgBa and AIB2 the 
peak popition is according to the inelastic neutron scatter- 
ing datatS. Open circles show the position ofrrtje-ifirst peak 
in the PC spectra for the corresponding metalstZHij. Straight 
dashed lines are guide to eye. 



Fig. H presents the recovered PC EPI functions (top 
panel) for the investigated diborides ZrB2 and NbB2 
calculated by Eq.(||) along with the High Resolutiaa. 
Electron-Encrgy-Loss Spectroscopy (HREELS) dataE£l 
for surface phonon dispersion (bottom panel). Details 
of the reconstruction of the PG-EPI functions according 
to Eq.(||) are described in Ref.Ej. 

Because of the large mass difference between the tran- 
sition metal and the boron atoms and strong covalent B- 
B bonds the boron modes are expected at much higher 
energy. For ZrB2 two additional maxima at 70 and 
100 mV are well resolved, while for NbB2 the high energy 
part of the spectrum presents a broad maximum around 



TABLE II: The phonon maxima and the EPI constant A in 
MeB2 compounds measured by PC spectroscopy. The fifth 
column shows thjB-ipipiximal energy for phonon features in 
the PC spectrumEZHij for the corresponding transition met- 
als: Me=Zr, Nb, Ta. 



Samples 


1"' peak 2 


peak 3^* peak Tiui^^^ 


Apc 




meV 


meV 


meV 


meV 




MgB2 


30(?) 


~ 60 


~ 90 


30 




ZrB2 


30 ±0.5 


68 ± 1 


94 ±2 


25 


0.06 


NbB2 


28 ± 2 


60 ±5 




28 


0.08 


TaB2 


20 ± 1 


40(?) 




20 


0.025 
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FIG. 2: The PC EPI function for ZrB2 and NbB2 recovered 
from the spectra in Fig.|l| (top paa^l) as compared to the 
HREELS data for these compoundgij (bottom panel). 



FIG. 3: Averaged for both polarities spectrum V2 oc cfV/dl^ 
after subtracting of a linear background of MgB2-Ag PC 
(Ro = 35 0, T=4.2K, Vite^=2.52 mV) in comparison to 
the neutron PDOS of MgB2t3. Inset: differential resistance 
K=dV/dl(y) of the same PC. 



60 mV. For TaB2 the high energy phonon peaks were dit. 
ficult to resolve, where according to a rough estimationEll 
the boron in-plane and out of plane displacement modes 
should have energies of 98 and 85 meV, respectively. No 
spectral features were found for the mentioned com- 
pounds above 100 meV. Q«r spectra are in line with the 
measured surface phononal3 for ZrB2 and NbB2 although 
the surface phonons, in general, use softer. For both com- 
pounds phonon dispersion studyllj (Fig.||, bottom panel) 
demonstrates a gap between 30-50 meV which separates 
acoustic and optic branches. Close to this energy region a 
minimum in our PC spectra occurs. Comparing the high 
energy parts of the ZijBe and NbB2 PC spectra, we may 
support the statement^, that for NbB2 the boron surface 
phonon modes are softer and more stretched than in the 
case of ZrB2 giving rise to structureless maxima around 
60-70 mV (see Fig.|]). Note that the upper boundaries of 
the ZrB2 and NbB2 PC spectra are at about 110 mV and 
90 mV (see Fig.g), what is much larger than their Debye 
temperatures of 280 K and 460 K, respectively, estimated 
from the Bloch-Griineisen temperature dependence of the 
resistivit}^. 

With the use of the EPI fimction the parameter A — 
2 ^ a^F{Lo)uj~^djjj wa^calculated. For many supercon- 
ductors it was foundllS that Apc — ^EUashberg- Table 
2 shows that Apc is rather low for the investigated di- 
borides what is in line with the small A values reported 
for PC-studies of the transition metal silicides NbSi2 and 
TaSi2E3 (Apc ~0.02). The possible reasons for a small 
A value are discussed in Ref.EHl. Our results show that 
NbB2 has the largest Apc among the studied compounds, 
therefore existence of superconductivity in NbB2 with ob- 
servable Tc is more likely. 



B. MgBa 

In Fig. I the V2 (eV) cx (fV/dP (eV) and dV/dI(y) 
(in the inset) characteristics are shown in the supercon- 
ducting state at zero field. In the energy-gap region the 
two gap features are clearly seen with Ai = 2.4 and 
A2 = 5.7 meV (see inset). Further on, we simply de- 
note by A the position of a dV/dl minimum. In Fig.^ 
the V2(V) oc d'^V/dP dependence after subtracting of 
linear backgrouniin comparison to the phonon density 
of states (PD0S)E3 is shown. Some accordance is seen 
between the positions of PDOS peaks and maxima in 
V2{V) dependence. 

The reproducibility of this kind of PC spectrum for 
different contacts is seen in Fig.^(b). Here, the contact 
resistance varies from 43 to 111 fl with gap minima at 
A = 2.1, 2.6, and 4.9 meV, respectively. The larger gap 
for curve 3 equals 7.0 meV. All the PC spectra (b) corre- 
late with PDOS. The slight variation is probably due to 
the anisotropy and different scattering rates in the con- 
tacts. Compared with the theoretical EPI functian-|£pr 
an isotropic one-band model (see, for example, RefJj'ffl), 
the peak at eV=60-;-70 meV of E2g boron mode is not 
too much higher in intensity, in accord with our earlier 
observationES. For the two-band model in clean limit, we 
expect an information mostly from the 3D band keeping 
in mind that we measure the c-axis oriented films. In this 
case, the theory predicts, that the EPI spectral function 
has the i?2g-peak intensity of the same order of magni- 
tude as the other peaks of PDOS (see Fig. 1 in RefBJ) 

Accordiag to our classification of the energy gap 
structureEj, it is due to the random orientation of the 
contact axis and scattering of charge carriers between 
two-band of Fermi surface, having gaps Ai and A2 
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FIG. 4: (a) Normalized to the normal state (at V ~ 30 mV) 
differential resistance at T=4.2K for 3 PC's with J?o = lll, 
43, and 45 f2 (curves 1-3 respectively). Curves 1,2 are shifted 
for clarity, (b) Second harmonic signal V2 with modulation 
voltages Vi(0)=2.5, 2.78, and 3.31 mV, respectively for same 
PC's. 
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FIG. 5: Normal state spectra of MgB2-Cu PC's taken at 41 K 
with Vi(0)=2.2mV. For spectrum No.l vertical arrows mark 
the low frequency bumps at about 20 mV. The horizontal dou- 
ble arrow bar stands for thermal smearing, determining the 
spectral resolution SV. 



mixed. With an increase of interband scattering the mag- 
nitude of Ai (A2) moves to the higher (lower) value, 
respectively. For dirty contacts, where the admixture 
of the 2D band is essential, only one gap maximum m- 
mains with a broad distribution around ~ 3.5 meVtj. 
In this case I is smaller than d (where / and d are the 
electron mean free path and the size of the contact, re- 
spectively), and the inelastic backscattering (jMitribution 
to the phonon structure, proportional to Z/ck£I, is small. 
Therefore, in the superconducting state, the observed 
phonon structure presents., mainly the elastic contribu- 
tion to the excess curren1n3. The elastic term is propor- 
tional to the energy dependent part of the excess current 
lexcieV), similar to the phonon struct ur e in the quasi- 
particle DOS for tunneling spectroscopyE3E3: 



dV ^ ' Ro 



A„, (eV) 



eV+J{eVr-AlieV) 



(3) 



where Ai„ (eV) is the gap parameter in the 3D band, 
induced by the 2D-band EPI. As seen in Fig. I (a), if the 
A value increases (curves 1-3), then the intensity of the 
phonon structure in the units of d{\n R) / dV = 2^/2V2/Vi 
increases too (Fig.^(b)). That is because the modulation 
voltage Vi decreases, whereas the amplitudes in V2{eV)- 
units remain approximately the same. Note, that for 
curve 3 in Fig. |(a) one has to take the lower gap Ai, 
because the current is mainly determined by 3D band, 
as mentioned above. 

In Fig.^, the normal state spectrum above Tc is dis- 
played together with the differential resistances in the 
superconducting state (see inset), showing the energy 



gap structure. In spite of an increase of the tempera- 
ture smearing above 20 mV at 41 K the residual phonon 
structure is still visible (curve 1 in Fig.||). The smeared 
phonon features are superimposed on the rising linear 
background. For curve 2, we see an increase in scattering 
at ~ 35 mV, where the acoustic phonon peak occurs, and 
the saturation around 100 mV, where the phonon spec- 
trum ends. In the normal state, only those nonlinearities 
remain, which are due to the inelastic processesEj. We 
stress that judging from the larger superconducting gap 
value (~ 3.5 meV for curve 2 in the inset) the essen- 
tial contribution is expected from the A2 {E) of 2D band 
for normal-state spectrum No. 2. This spectrum should 
contrast with curve 6 in Fig.|^ (see below), where the di- 
rect contribution from the 2D band is small, due to lower 
value of A. Thus, the backscattering processes from the 
3D band are mostly essential for spectrum 6 and the 
phonon features are not resolved at high energies. The 
shape of spectrum 1 in Fig. ^presents the behavior inter- 
mediate between these two extremes, although its energy 
gap is approximately equal to curve 6 in Fig.^. Beside 
the visible phonon features in the range 30 100 mV, it 
possesses a low energy bump at about 20 mV. 

The phonon spectra of PC with a small value of energy 
gaps are characterized by the presence of low frequency 
phonon peaks. The small peak at energy of about 25 mV 
(Fig. ^ cuKiies 1, 2) is visible, where a tiny knee exists on 
the PDOSta (see PDOS in Fig.|). In the normal state 
(at T > Tc), these low frequency peaks transform into 
the S-shape structure in d^V/dP{V) spectra (curve 6), 
corresponding to the wide minimum of dV/dI{V) near 
zero bias. This low-frequency structure is hardly due 
to the remnants of superconducting quasigap at T > Tc, 
since it is absent injunction No. 2, whose characteristic is 
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FIG. 6: Second harmonic dependences on field and temper- 
ature averaged for botfi voltage polarities. The temperature 
and magnetic field are: 4.2 K, OT; 4.2 K, 4T; 10 K, 4T; 20 K, 
4T; 30 K, 4T; 40 K, 4T, for curves 1-6, respectively. The 
modulation voltage Vi(0) is 2.2 mV. A=2.7meV. Curves are 
shifted vertically for clarity. 

shown in Fig.||. Rather, it could be thought of as strong 
interaction in the 3D band with very low-frequency exci- 
tations, whose origin is not clear yet. On increasing the 
temperature, the low frequency peak broadens further 
like common spectral features do in the normal state. 

Fig. H shows (compare curves 1 and 2) that the field 
smears out the intensity of the high energy peaks, which 
are induced in the 3D band by EPI of the 2D band. 
The disappearance of phonon peaks at field and tem- 
perature rise proves that they do not belong to the in- 
elastic backscattering processes, which should have the 
same intensity both in the superconducting and in nor- 
mal states. It seems more plausible that the high energy 
phonon peaks are due to the elastic contribution in the 



excess current (^) induced by EPI from the 2D band, as 
was already stated above. 

Conclusion. We have measured the PC spectra in tran- 
sition metal diborides: ZrB2, NbB2 and TaB2. For all 
compounds the main phonon peak position was estab- 
lished, the PC EPI function was recovered, and the EPI 
parameter A was determined. The obtained small A val- 
ues question strongly the reported bulk superconductiv- 
ity in these compounds with remarkable Tc. To draw a 
more weighty conclusion about details of EPI and the A 
values in the presented MeB2 family, a theoretical calcu- 
lation of a|Q-F(a;) with the mentioned A'-factor and their 
comparison with experimental data is very desirable. 

The PC EPI spectra of the superconducting MgB2 dif- 
fer even qualitatively from those measured for the men- 
tioned transition-metals diborides. For MgB2 the repro- 
ducible peaks on <PV/dI'^{V) are seen in the supercon- 
ducting state, which disappear in the normal state. It 
proves that they are due to the energy dependence of the 
superconducting order parameter. The superconductiv- 
ity here is presumably due to the 2D-band EPI, which 
induces the A(i?) structure in the 3D band. At a small 
value of the superconducting gap the phonon structure 
is weak, and its intensity begins to increase with growing 
A. The robustness against the magnetic field also grows 
notably with increasing A. 

In the normal state the intensity of the inelastic spec- 
trum also correlates with the value of the superconduct- 
ing gap. The tendency is observed that the high energy 
phonon contribution to the PC spectra becomes more 
prominent for a larger gap. This is supported also by 
our recent data on MgB2 single crystalsE^. For a smaller 
gap, the low- frequency bump appears in the spectra. It 
might be due to EPI with some unknown low-frequency 
excitationE3, but further work is needed to exclude that 
it is not caused simply by phonon peaks of the normal- 
metal counter electrode. 
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